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In the commented paper [PRD 78, 034032 (2008); larXiv:0803.2956] Deandrea et al. claim that 
they provide a model- independent determination of the quantity Bo(md — m u ) based on the analysis 
of the decay of eta to three neutral pions that takes the strong-interaction correction as an unknown 
parameter. In the following we show that the two constraints presented there as independent 
equations are in fact two approximations of a single constraint connecting together Bo(md—m u ) with 
the strong-interaction correction to the r\ —¥ 3tv° decay amplitude. Thus, without some additional 
information on the latter (such as its correct value at least at one kinematic point) it does not 
lead to any determination of Bo(rrid — m u ) and the numerical values presented there should not be 
taken seriously. We have also found that the numerical result for the studied quantity obtained in 
the commented paper is just an accidental number stemming from some numerical errors that had 
occurred in the computations presented there. 

PACS numbers: 12.39. Fe, 11.30.Rd, 14.65. Bt, 13.25. Jx 



I. INTRODUCTION 

The Bose symmetry of the decay amplitude of the 
rj — > 3tt° decay translates into a very simple form of the 
amplitude that is valid up to the two-loop ordeiQ 

M(s, t, u) = M(s) + M{t) + M{u). (1) 

The Mandelstam variables s, t, and u fulfill the constraint 

s + t + u = m„ + 3m 2 = 3s . (2) 

This decay violates the isospin symmetry. There are 
two sources of the isospin breaking, the m^ — m u mass dif- 
ference and the difference of the electric charges of these 
two quarks. The considered decay amplitude thus has 
two types of contributions, the strong one proportional 
to rrid — m u , and the electromagnetic contribution pro- 
portional to e 2 . Neglecting the contributions of orders 
(rrid — m u ) 2 , e 4 [and all higher ones], which should be 
suppressed with respect to those of the lower orders, we 
can write the s-wave amplitude in the forir0 

M(s) = -TT-J^J (1 + M«) + <W0O) + ^em(s). (3) 

This parametrization employs the low-energy constants 
of the lowest order of chiral perturbation theory, the pion 
decay constant F n and the scalar quark condensate Bq 
hidden in 

X=B Q (m d ~m u ). (4) 
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1 For a more complete discussion of this amplitude, we refer to [l| 
and 0. 

2 Note that this form can be kept valid even if we include the 
neglected higher orders by redefinition of the meaning of the 
individual Ss. 



The naive power-counting X^str(s) = 0(p 2 (md — ra u )), 
X<5cm(s) = 0(e 2 (m d - m u )), 5 em (s) = 0(p 2 e 2 ), together 
with the explicit computation in the first orders in chiral 
perturbation theory leads us to expect that the electro- 
magnetic contributions are tiny with respect to the strong 
one, namely 

8 e m(s) < X<W(s) and S em (s) < S siI (s). (5) 

In addition, the strong contributions of higher chiral or- 
ders are expected to be significantly smaller than the 
leading order contribution, i.e. 

Sstr « 1. (6) 

This decomposition of the decay amplitude shows how 
its computation can be useful for a determination of the 
quantity \ — ^o(™d ~ m u) giving us information about 
the individual m u and quark masses. Indeed, pro- 
vided we know all of the Ss appearing in ([3]), from the 
total decay rate of the considered process we can deter- 
mine the unknown \. 

An additional independent determination of x comes 
from the comparison of this amplitude squared with its 
experimentally measured energy dependence. In the ex- 
isting experimental studies this is performed by a deter- 
mination of the parameters of its expansion around the 
center of the Dalitz plot s = t = u = sq, the so-called 
Dalitz plot parametrization, 

\M(s,t, u)\ 2 - |X(.s )| 2 (l + 2az + ...), (7) 

where z is a dimensionless parameter denoting the dis- 
tance from the center of the Dalitz plot 

z = ^((s-so) 2 + (t- S o) 2 + (u- So ) 2 ) (8) 

with 

Q = mi(m v - 3m ff o ) 2 . (9) 



2 



Note that the latter determination is possible only if we 
fix the physical normalization of the amplitude by 5 em (s). 
If we neglect this contribution, the desired quantity \ 
appears in the overall normalization of the amplitude ([3]) 
and even a precise measurement of the round bracket in 
(|7|l would not enable a determination of x- 



this expansion, we compute the absolute value squared of 
this amplitude. By employing the form (J3J) of its s-wave 
part and by comparison of the result with the Dalitz plot 
parametrization ([7]), we obtain "the first relation" 

X 2 («str + "em) + X^cmA^ = X 2 a, (12) 



II. KEEPING THE STRONG-INTERACTION 
CORRECTION AS AN UNKNOWN 
PARAMETER 

Having these two independent methods of the determi- 
nation of x, there arises the question whether it is possi- 
ble to combine them in order to determine this quantity 
without any knowledge on the 5 s (s). 

If S s (s) was a constant, we would indeed have two inde- 
pendent equations for two unknowns. However, <5 s (s) is a 
nontrivial function of s. In each of the methods there ap- 
pears a different characteristic of this function — in the 
determination that uses the decay rate of i] — > 3n° we 
have the integration of S s (s) over the physical kinematic 
region, whereas in the method using the Dalitz plot pa- 
rameter a, the expansion of 5 s (s) around the point s = sq 
is used. Even if we neglect all the higher derivatives of 
<5 s (s), in addition to searching for a correct value of the 
constant S s (sq), we would need to solve the equations also 
for the second derivative 5"(so), whose real part is in [TJ 
hidden in a s . 

We can therefore expect that the equations presented 
in [l| as independent equations in which there appear just 
a s and B^md — m u ) are in fact dependent and thereby 
insufficient for the determination of both of these un- 
knowns. This will be shown in the rest of this comment. 



where we have denoted the real parts of the second 
derivatives of <5s according to 



a st Y = — Re<5" tr (s ), a c 



■ReO* ), ( 13 ) 



(ScmA^ = —F% -j= Re S" m (so) 



(14) 



and neglected all terms of orders 0(x 2 S 2 ), 0(S^ m ), 
0(xS§em), 0(x 2 Sa), and 0(x5emOt)- [In this counting of 
orders, under a shortcut 0(S) we understand both differ- 
ent orders 0(S stI ) and 0(5 cm ) together. We also assume 
that each S"(sq) is of the same order as <5(so).] 

If we assume the validity of ([5]) and © and use the 
smallness of the experimental value of a, it is reasonable 
to neglect these orders and with this relation we have 
come to the same conclusions as the authors of [l[. The 
only difference is in the numerical evaluation of a om as 
will be discussed below. 



B. The relation that uses the decay width 

The decay rate of rj — > 3tt° can be obtained by the 
integration of the amplitude AA{s,t,u) over the whole 
decay kinematic region, 



A. The relation following from the Dalitz plot 
parametrization 

The first relation follows from expanding the decay am- 
plitude M(s, i, u) around the center of the Dalitz plot by 
using ([1]) and 



M(s,t,u) = — 



3=0 



(— •<»)£ + (*- «o)| 



+ (u - s ) 



d_ 

du 



M(s,t,u) 



(10) 



Since -j^M(s,t,u) = J^M(s), the symmetry of the am- 
plitude leads to vanishing of the term with the first 
derivative [because of ©] and the term with the second 
derivative contains z from (|SJ), 

M(s,t,u) = 3M(s ) + QzM"(s ) + 0(z 3/2 ). (11) 

With 0(z 3 / 2 ) we formally denote the terms of the higher 
than second order in (s — sq), (t — sq) or (u — sq). From 



r = i 1 



ds / dt \M(s,t,u)\ 2 , 



6 2567r 3 m3 J 4m 2 o j t (s) 



where the boundaries of the t-integration are 



(15) 



t± (s) = \ (3s - s ± X^l (sKo (s)) (16) 

and A and a are the usual Kallen triangle function and 
the kinematic square root, respectively (cf. 0). 

Note the appearance of the combinatorial factor 1/3! = 
1/6 connected with the presence of three indistinguish- 
able particles in the final state. This factor was omitted 
in [l|. In order to reproduce the results of [l| as closely 
as possible, we include this factor into a redefinition of 
r, i.e. from now on with T we denote the experimental 
value of the decay rate multiplied by the factor of 6. 

We can exchange the order of integrations in ([TSl and 
thanks to the symmetries of the kinematic region the 
boundaries of the internal s-integration would then be 
t±(t) depending on the variable t of the external inte- 
gration. Therefore, by interchanging the name of the 
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integration variables, we obtain 

(m v -m !r o) 2 /■*+(») 

ds / dtA(s,t) 

4m 2 „ Jt-(s) 

ds / dtA(t,s) (17) 

4m 2 •'t-W 

for any function A(s,t). Moreover, we can perform the 
same changes also with the variable u and thus the inte- 
gration of any function A(s, t, u) gives the same result as 
the integration of this function with any of the Mandel- 
stam variables interchanged. 

By using this fact, the decomposition ((3|), and neglect- 
ing again the orders 0(x 2 S 2 ), 0(S 2 m ) and 0(xSS em ), we 
obtain the decay rate in the form, 



T = X 2 (7trco + 7str + 7cm) + X7cm, 



(18) 



where 7s are introduced in accordance to [1] as the fol- 
lowing integrals 

2 2 
7str = j^4-F[Re 4tr(s)] , 7cm = ^4 J" [Re <5 om (s)] , 



7tr 



The functional T denotes the integration 



(19) 

J-[Re5 em (s)]. (20) 



1 



2^,vM/3./ 4M2 dS / dtfis) 



1 



256tt 3 M3 



(il^-M^) 2 



dsXl%)a^(s)f(s). (21) 



Now, we introduce a symbol T for the decay rate of 
this process in the situation the amplitude was constant 
and everywhere equal to its physical value at the center 
of the Dalitz plot, T — |7W («o ) 1 2 ~^"[1] - I n order to be more 
suggestive, we divide it into three parts 



X 2 (it 



Tstr Torn) XTc 



(22) 



with 

2 2 

7str = ^lRe4tr(so)^"[l], 7cm = oT77 Re <S G m (so )-^[l] , 



3Fi 



3Fi 



(23) 



7tr 



7tree, 7c 



V3^ 2 



Re<J cm (a )^[l]. (24) 



With the help of this quantity, the decay rate can be 
written as 

r = r + ^ T [Re (<5 str (s) - Wo))] 



2 X 2 



3F4 
(>X 



J 7 [Re (<5 em (s) - <5cm(s ))] 



V3F2 



Re 



(&m(«) - 4m(s ))] , (25) 



where we have neglected again the same orders as in (|18l) . 

On the other hand, the integration of the Dalitz plot 
parametrization ([7]) leads to 

T = |-M(s )| 2 (V[l] + ^aF[(s - s ) 2 ] + F[0(z 3 ' 2 ) 
= r + ^- a \M(s a )\ 2 T[(s - s ) 2 ] + T[0(z 3 ^ 2 )]. (26) 



If we set the higher Dalitz plot parameters included in 
the dots of ([7]) equal to zero or if we find JF[0(z 3 / 2 )] negli- 
gible by any other reasoning, the comparison of equations 
(f25"l) and (|26|) gives the relation 



X 2 «str + X 2a em 



XOem(s) = ~a e 



where we have introduced the symbols 

2Q T[Re{6 stT (s)~6 stv {s ))} 
9 .F[(s-so) 2 ] 
2Q T [Re (S em (s) - <5 cm (s ))] 



Qstr 



(27) 

(28) 
(29) 

(30) 
(31) 



From the first line of (|26|) . we can express the value of 
the amplitude at the center of the Dalitz plot and obtain 



J-[(s- S0 ) 2 ] 



^QF 2 



Re 



V3 FKs-so) 2 } 
-3F*a\M(s Q )\ 2 . 



QT 



+ 9aT[(s - s ) 2 



(32) 



which is exactly its definition of [l|. Up to the identi- 
fication of a str = a s tr performed in [l|, in (j2"?| we have 
obtained the "second independent equation" for x- 

We should recall that this relation is an approximate 
one. Similarly, as in the first relation we have neglected 
the terms of orders <3(x 2 <5 2 ), 0(6 2 m ), 0(xS5 cm ). But 
in this relation we have also made an assumption about 
the higher order terms of the Dalitz plot parametrization. 
Setting them to zero is equivalent to setting the third and 
all the higher derivatives of the amplitude equal to zero. 
We can also understand this assumption as neglecting of 
the terms stemming from such higher derivatives. 

This can be done also for the individual contributions 
from the decomposition ([3]), i.e. we can expand all of the 
5s around the center of the Dalitz plot, 

5( s )=5(s )+(s-s )5'(s )+^(s-s ) 2 6"(s )+... , (33) 

and neglect all the terms appearing within the dots. 

Thanks to the symmetries of the integration over the 
decay kinematic region [discussed around (1171) ]. the first 
derivative does not contribute, 



F[s - s ] = -T [(s - s ) + (t - s ) + (u - s )] = 0, 



(34) 
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and from expansion (|33l) there follows 

T [6(a) - S(s )} = ±8"(8 )F[(a - s ) 2 } + .... (35) 

In other words, if we neglect these higher-order-deri- 
vative terms, we get 

As we have discussed above, in [l[ the first of these equa- 
tions was used. However, if we neglect the third-order 
derivative of <5 s tr(s), there is no reason for keeping the 
third-order derivatives of <5 em (s) and 6 em (s), as these 
functions should be considerably smaller than S s t r {s)- 

With that we have identified the left-hand sides of both 
of the relations (fT2]) and (|27|) , so we already see that they 
are not independent. Now, let us proceed to the right- 
hand sides. We start with the definition of a exp (|3ip and 
instead of evaluating |.M(so)| 2 from the experimental T 
through (|26|) . we use once again the decomposition ([3]) 
[which is in (j2"T|) used anyway] . By that we obtain^ 

a cxp = -X 2 " + 0{ax 2 S) + 0(ax<5 cm )- (37) 

After using this identification, relation (1271) goes exactly 
into relation (fT2|) . 

In conclusion, each of these relations is a different ap- 
proximation of the same unique relation between the 
function 6 s t T (s) and the quantity the first one using 
the smallness of a, whereas the second one is neglecting 
higher-order derivatives of the amplitudes. These differ- 
ent omissions cannot lead to two independent relations 
for the unknown parameters since their use means that 
we assume the reasonableness of both of the omissions, 
which should be thus used together; and this leads to an 
unique approximate relation between the parameters. 



III. CONCLUSIONS 

We have shown that the "independent parameters" 
a cm , a cm and a om ,a em should fulfill the approximate 
relations (131)1) (up to the contributions of higher-order 
derivatives, which should be numerically small). Indeed, 
if we evaluate the expressions for 5s from we obtain 
an agreement between the values. 



We should therefore understand these relations in the 
following sense. We have a constraint between \ an( i ai l 
the functions 6(s), which is in the form of (|27| . where we 
have two different ways how to approximate the param- 
eters di from its left-hand side, either through (|28l) - (|30|) 
or through (|36|) . From these two approximations, we 
can determine the values of a, and estimate their errors 
stemming from the higher-order corrections neglected in 
these two different enumerations. Similarly, the right- 
hand side of this relation, the parameter a oxp , can be ob- 
tainecfl either from (j3"2")) or from ([57)1 . In order to solve 
this constraint and obtain the quantity B()(md—m u ) from 
it, we need some additional assumption or constraint on 
the physical amplitude M.(s,t,u) or on its part <5 s tr(s)- 
An example of such an additional condition is to set the 
value of <5 s tr(s) at some specific point, where we for in- 
stance believe its chiral computation (cf. Q). 

Why have the authors of [l| not found these conclu- 
sions in their numerical analysis? We have checked the 
results of [l| and found an error in the expression given 
there for aim, which does give the (wrong) numerical 
value of this quantity that has been presented in 
Similarly, the numerical values of a, a eX p and the value 
of x should approximately [cf. footnote 4] fulfill relation 
(l37|) . which is not the case for the values presented in 
As we have discussed above, the reason is that in 
the relation ([52)1 one should introduce not the experi- 
mental value of r but this value multiplied by 6. Both 
of these errors shift the coefficients of relation (1271) by 
a quantity of the order of the correct coefficients, and 
such a shift comprises an unphysical second independent 
constraint on a str and \ in addition to the unique inde- 
pendent approximate relation (|27[) . This explains quite a 
reasonable magnitude of the values for \ = Bo(m,d — m u ) 
and Q s tr obtained from such an analysis in [l[ . We stress 
once again that these numbers are chosen from the whole 
two-parametrical space of solutions of relation (f2"Tj) just 
by accident and one should not take them seriously. 
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3 Note that thanks to the smallness of a and expecting the validity 
of (O and JB), we can use here [similarly as in relation JT2}] just 
the leading order of M (s) with all <5s from decomposition J3} 
neglected. Exactly this is how we have got rid of the unknown 
5str(«o) — however, by that we have also got rid of the genuine 
second independent relation [dependent on this <5 s tr(so)] that was 
discussed in the introduction. 

4 According to the author of this comment, the determination (132 I t 
(with the correct definition of T including the factor of 6) is 
more reasonable as it ignores just an integral of higher Dalitz 



plot parameters, which should be small due to the smallness 
of the decay kinematic region of this process. In fact, after a 
measurement of the higher Dalitz parameters, one can simply 
make its value more precise by adding these parameters into 
H26I I. On the other hand, the determination H37II rests heavily 
on the condition that 5s in J3|l are very small. However, the 
comparison of these two determination by using the value of Bo 
from the recent lattice determinations and of — m u from Q| 
leads to ^fz^i ~ 1-9, i-e. <5s can be significantly larger than 
expected and one should be careful when employing J5{ and (5) . 
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